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Density Functional Theory for Superconductors 
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We apply the density functional theory for superconductors (SCDFT) based on the local-density 
approximation (LDA) to alkali-doped fullerides A^Cqo with the face-centered cubic structure. We 
evaluate the superconducting transition temperature (T c ) from first principles considering energy 
dependence of electron-phonon coupling, the mass renormalization, and the retardation effect. The 
calculated T c =7.5, 9.0 and 15.7 K for A=K, Rb, Cs are approximately 60 % smaller than the 
experimentally observed values. Our results strongly suggest necessity to go beyond the framework 
of the Migdal-Eliashberg theory based on the LDA. 

PACS numbers: 



-Introduction. The doped fulleride superconductors 
A 3 C 60 (A =alkali metal) [J Q, which exhibit maxi- 
mum transition temperature (T c ) of 40 K, have pro- 
vided a fertile playground for theoretical and experimen- 
tal studies. The most significant feature of the fullerides 
is the narrow metallic bands formed by molecular or- 
bitals, whose energy scale compete with the vibrational 
frequencies and electron-electron interactions. Moreover, 
recent experiments revealed that the T c -volume (V) curve 
for this series shows a dome-like dependence near the 
superconductor-Mott insulator transition [H-Q . This de- 
pendence is, similarly to the celebrated superconducting 
dome in cuprates, reminiscent of a crossover from weak 
to strong correlation in this system. 

Motivated by these properties, various theoretical 
studies have investigated unconventional pairing mech- 
anisms [l], [|| . On the other hand, there have also been 
a received idea that the superconductivity in this system 
is explained by the conventional phonon-mediated pair- 
ing mechanism. Full s-wave gap with spin-singl et p air- 
ing C-isotopc effect coefficient of >0.20 [Hf[ 
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and the coherence peak in the NMR and ^(SR relaxation 
rate 15, 16[ have been experimentally observed. In par- 
ticular, in the T c -V plot, the regime where T c and V 
positively correlate is seemingly consistent with the BCS 
theory; increasing V results in smaller bandwidths, larger 
density of states (DOS) at the Fermi level, and subse- 
quently stronger electron-phonon coupling. Hence, the 
applicability of the phonon mechanism is still unsettled. 

The Migdal-Eliashberg (ME) theory [It}, which is a 
quite widely applicable theory of phonon-mediated su- 
perconductivity, includes the self energy with the lowest- 
order exchange contribution of the dressed phonons and 
the static screened Coulomb interaction. The ME theory 
described by the Kohn-Sham orbital based on the local- 
density approximation 0, [l9[ (KS-LDA) has enabled us 
to consider fine details of materials (20I - I22I ] . Moreover, 
the recently developed density functional theory for su- 
perconductors (SCDFT) [23, 24 1 has provided us a way to 



calculate T c based on the ME theory nonempirically. The 
SCDFT treats effects of the interactions such as the mass- 



renormalization 17 and the retardation effect |25| taking 
the detail of the electronic structure. The T c calculated 
with the SCDFT has shown remarkab ly g ood agreement 

However, 



24. 26. 2 



with experimentally observed T c 
its applications to molecular solids has not been reported 
due to its expensive computational cost. In this paper, 
we apply the SCDFT to fee A 3 C 60 [A=K and Rb under 
ambient pressure (T c =19 and 29 K), and Cs under the 
optimum pressure of 7 kbar (T c =35 K)] focusing on the 
regime where T c and V positively correlate. We calculate 
T c to see if the SCDFT reproduces the absolute values 
and the alkali-metal dependence of the experimentally 
observed T c , with which we examine the applicability of 
the ME theory with the KS-LDA in the present system. 
The calculated T c suggests that we need to consider some 
factors missing in the framework of the ME theory based 
on the KS-LDA. 



-Method. In the current SCDFT [23j, l2j we solve the 
gap equation given by 
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Here, n and k denote the band index and crystal mo- 
mentum, respectively, A is the gap function, and f3 is 
the inverse temperature. The energy i? n k is defined as 
i?rik=v / C 2 ik + an d £nk = £nk — /it is the one-electron 
energy measured from the chemical potential fi, where 
e n k is obtained by solving the normal Kohn-Sham equa- 
tion %Ks|^«k) = EnklVnk) with Hks and |<£ nk ) being 
the Kohn-Sham Hamiltonian and the Bloch state, re- 
spectively. The functions Z and /C are the exchange- 
correlation kernels describing the effects of the inter- 
actions. The kernels describing the standard electron- 
phonon mechanism, K,=JC ph +IC el and Z=Z ph , have been 
proposed 



23l . l24j . Namely, the phonon contributions 



(K. ph and Z ph ) were formulated referring to the ME the- 
ory, and the electron contribution (/C cl ) corresponds to 
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the screened static Coulomb interaction scattering the 
Cooper pairs. 

Since the fullcridc superconductors involve high- 
frequency phonons, the electron-phonon interaction has 
strong dependence on both £„k and £ n 'k'- In order to 
treat this effect, we use the nk-resolved form for IC ph and 



Z? h defined by Eqs. (9) and (11) in Ref. |2J, which re- 
quire the electron-phonon matrix elements , k , and 

the phonon frequencies w„ q as inputs. For /C cl , on the 
other hand, we use the form given by Eq. (13) in Ref. I29I 
which is based on the static random-phase approxima- 
tion (RPA) [28[ and properly treats the local-field effect 
due to the spatial dependence of electron density. 

-Result and discussion. We calculated the band struc- 
ture, phonon frequencies, electron-phonon and electron- 
electron interactions, and T c for fee A^Cqq with A=K, 
Rb and Cs. All of our calculations were performed 
within the local-density approximation using ab initio 
plane-wave pseudopotential calculation codes Quantum 
Espresso (30l - [33 | . Input structural parameters [36| were 
determined b y fu ll relaxation and the Murnaghan equa- 
tion of state [37|. For A=Cs, we considered the case of 
the optimum pressure of 7 kbar Q. Phonon frequency 
and electron-phonon coupling [38[ were calculated by the 
density functional perturbation theory [3j|. The dielec- 
tric function used for /C el was calculated within the static 
RPA. The SCDFT gap equation [Eq. was solved 
with the random sampling scheme given in Ref. l4oL with 
which the sampling error in the calculated T c was ap- 
proximately 3%. 

Let us start from the calculated DOS of the partially 
occupied t\ u bands in Fig.[TJ The general trend is consis- 
tent with the previous calculation based on the general- 
ized gradient approximation and the experimental lattice 
constants [13] • As anticipated previously 0, H| , replac- 
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FIG. 1: (Color online) DOS around the Fermi level. Inset is 
the view in a broader energy scale, where the characters of 
the bands are specified. 



ing lighter alkali metal elements with heavier ones (from 
K, Rb to Cs) leads to slightly larger DOS at the Fermi 
level AT (0) (see Table fl]). More significantly, we also see 
that the bandwidth becomes narrower. The relation be- 
tween these changes and the electron-phonon coupling is 
discussed later. 

Table [TT] summarizes our calculated frequencies of 
the L-point ffg-derived modes, which are distinguished 
as five-fold degenerate branches with strong electron- 
phonon coupling. The experimentally observed and pre- 
ceding theoretical frequencies are also given for compar- 
ison. The agreement between our calculation and ex- 
periments is extremely good, which illustrates that our 
calculation properly describes the phonon properties of 
the present system. The alkali-metal dependence of the 
frequencies is little, which is due to the intramolecular 
property of the modes. 

We next show the T c calculated by the SCDFT with 
only the phonon contributions to the gap-equation ker- 
nels (K. ph and Z ph ). The calculated T c (red solid 
square) is higher than the experimental T c , which 
is because of the absence of the electron contribu- 
tion. These values are consistent with the recent 
calculation based on the Eliashberg equation [l7| by 
Koretsune and Saito [471 ] . Interestingly, the exper- 
imentally observed alkali-metal dependence is repro- 
duced. In order to examine the origin of this depen- 
dence, we calculated the electron-phonon coupling co- 
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efficient A 



AT(0) 



JV(0) ^kqrm'i/ 



and the characteristic frequency u)y n 



£>(£nk)<5(£ri'k+q) 



n(o)= exp 



N(0)X N(0 ) 

„'k+q)lnw„ q j. The calculated 
values are listed in Table HI By replacing lighter alkali- 
metal elements with heavier ones, Ajv(o) is slightly en- 
hanced due to the increase of N(0). However, when we 
substitute Aat(o) and winjy(o) into the McMillan- Allen- 
Dynes (MAD) formula j48|, T c =ff exp[-1.04(l + A)/A] 
(with the Coulomb pseudopotential fi* set to 0), the de- 
pendence of the resulting T c (blue open circle) is not as 
significant as that obtained from the SCDFT. Alterna- 



tively, we calculated Xn(^) an d <*>\ n 



MS 



defined by the 



TABLE I: Calculated parameters representing the electronic 
structure and the electron-phonon and electron-electron in- 
teractions. 





K3C60 RbaCeo 


CS3C60 


A/(0)[/(eV spin) 


8.352 


8.609 


9.328 


^iV(O) 


0.562 


0.570 


0.603 




0.489 


0.542 


0.652 


<^ln,N(0) [K] 


1071 


1054 


1052 


w ln.iV(5) [K] 


932 


944 


940 


Z 


0.350 


0.367 


0.396 


A* 


0.379 


0.370 


0.362 
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TABLE II: Experimentally observed and theoretically calculated T-point phonon frequencies (cm 1 ). The dashes denote the 
splitting induced by the crystal field. 



expt. 






present 








theory 


a/— i 


K3C60 


K 3 


C60 


RbaCeo 


Cs 3 


C60 


T^Ceo K3C60 


H e (l) 273 


271 


262- 


-271 


261-269 


261- 


-270 


281 


252-258 


H g {2) 437 


431 


422- 


-422 


420-422 


418- 


-421 


454 


407-404 


flg(3) 710 


723 


685- 


-689 


686-688 


687- 


-689 


753 


658-663 


flg(4) 774 




779- 


-779 


779-780 


779- 


-783 


785 


737-740 


H g (5) 1099 




1111- 


-1116 


1111-1116 


1113- 


-1120 


1091 


1019-1023 


flg(6) 1250 




1268- 


-1274 


1268-1273 


1271- 


-1275 


1290 


1137-1136 


H g (7) 1428 


1408 


1403- 


-1408 


1402-1405 


1406- 


-1407 


1387 


1349-1348 


flg(8) 1575 


1547 


1532- 


-1537 


1532-1536 


1532- 


-1538 


1462 


1532-1530 



a Raman scattering measurement, Ref. 43 . 
6 Raman scattering measurement, Ref. 44 . 

c Ab initio LDA full-potential calculation based on linearized muffin-tin orbital method, Ref. 45. 
d Ab initio LDA pseudopotential calculation based on the mixed basis method, Ref. 46 . 



following formulae |49|, [5 
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These formulae explicitly treat the energy conservation 
in electron-phonon scattering, and therefore include the 
effects of the electronic states within the phonon energy 
scale; since the scattering involves energy exchanges of 
order <0.2 eV, electronic states within this energy range 
should contribute to the pair formation. As a result, the 
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FIG. 2: (Color online) Calculated T c s: Solid squares denote 
the values calculated using the SCDFT gap equation with 
only the phonon contribution (K. ph and Z ph ), and open (solid) 
circles denote the values derived from the MAD formula (see 
text) using X N(0) (Ajy(o) and w in,JV(o) (wi n ,jv(o) in Table U 
The triangles represent the experimentally observed values. 



dependence of the calculated Ajv(^) is more noticeable 
than that of Ajy(o) i an d the corresponding T c derived from 
the MAD formula (blue solid circle) well reproduces the 
dependence of the T c calculated by the SCDFT and the 
experimentally observed T c . The present analysis clar- 
ifies the significance of the electronic states within the 
finite energy range, not only at the Fermi level. 
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FIG. 3: (Color online) Calculated gap function for Cs 3 C60 
under pressure of 7 kbar with T=0.01 K. The characters of 
the three bands are specified. 



We also found an important 
renormalization factor Z=Z^\^ n] 
In usual cases, Z is as large as Ajv(o) [2J|: but our cal 



aspect of the mass- 
_>o given in Table Q] 



culated Z is much smaller than A 



N(0) 



or A 



This 



is because that the t\ u bands are energetically isolated 
from other bands. The main contribution to the mass 
renormalization around the Fermi level generally comes 
from electron scattering to the states distributed within 
the several times of the Debye frequency. In the present 
case, however, the energy scale of the Debye frequency is 
as large as the bandwidth of the ti u bands, and there is 
no scattering channel in the gapped region (see the in- 
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set of Fig. [T]). This weak mass renormalization results in 
relatively hi gher T r than expected from the conventional 
calculations jl7l.l48j. 

Next let us move on to the results obtained 
with the electron contribution (/C cl ). The strength 
of /C cl is represented by its Fermi-surface average 
M=ivIo)S«k«'k /C «kn' k '5(e i k)<5(^kO (Sec Table H}. We 
display in Fig. [3] the gap function in T=0.01K for A=Cs 
under pressure of 7 kbar. The values of the gap function 
in the ti u states are positive, whereas those in the highest 
doubly occupied h u and the lowest unoccupied ti g have 
negative sign. Such sign inversion of the gap function in 
the high energy region represents the retardation effect 
in the SCDFT [24[. Here, the absolute values in the high 



mately 30 % by increasing the exchange contribution in 
the sclf-consitent calculation of the wavefunctions 51 



energy region are quite comparable to those in the low- 
energy region, which signifies the strong retardation ef- 
fect c omp ared with the previously reported conventional 
cases [24], [26J, |27j . This is due to large intcrband electron- 
electron Coulomb interaction. 

Finally, we show the calculated T c for A=K, Rb, and 
Cs in Fig. |4] together with the experimentally observed 
T c . Thanks to the energy dependence of electron-phonon 
coupling, the alkali-metal dependence of the experimen- 
tally observed T c is well reproduced. In spite of the weak 
mass renormalization and the significant retardation ef- 
fect, the absolute values arc 7.5, 9.0 and 15.7 K, which 
are approximately 60% lower than the experimentally ob- 
served T c (19, 29 and 35 K). The present underestimation 
is in clear contrast with the previous applications to the 
conventional superconductors [II EE}. 

Since the phonon frequency scale in the present sys- 
tem is quite large, the theoretical T c is sensitive to the 
input interactions. In fact, with „/ k /| 2 multiplied by 
1.2 (0.8), we obtain T c =17.5 (1.5), 20.6 (2.3), and 31.7 
(5.0) K for A=K, Rb and Cs, whereas we obtain T c =5.8 
(9.6), 7.7 (11.3), 14.7 (18.0) with IC cl multiplied by 1.2 
(0.8). Concerning the ab initio calculation of the inter- 
actions, on the other hand, a recent paper reported that 
the electron-phonon interaction is enhanced by approxi- 
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FIG. 4: (Color online) Calculated T c by solving the SCDFT 
gap equation with the electron contribution /C cl compared 
with the experimentally observed values. 



These facts imply a possibility to fill the discrepancy be- 
tween our results and experiments by considering the fea- 
tures neglected in the present ME theory based on the 
KS-LDA. 

-Summary and conclusion. Using the SCDFT, we per- 
formed non-empirical calculations of T c in fee A3C60 
(A=K, Rb, Cs). We focused on the energy dependence 
of electron-phonon coupling, the weak mass renormaliza- 
tion and the strong retardation effect. Our calculated 
values of T c were 7.5, 9.0 and 15.7 K for A=K, Rb and 
Cs (under pressure of 7 kbar), which are approximately 
60% smaller than the experimentally observed values (19, 
29 and 35 K). The present results indicate a necessity to 
go beyond the ME theory based on the KS-LDA even for 
the regime where T c and V positively correlate. 
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